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The  Application  of  a  Statistical  Analysis 
Software  Package  to  Explosive  Testing 


1.  Introduction 

Many  ordnance  evaluations  involve  the  study  of  reliability  and  safe  functioning  of  an 
energetic  material  or  device.  This  calls  for  statistical  analysis  methods  to  accurately  evaluate 
a  given  set  of  data.  For  explosive  testing  purposes,  Material  Research  Laboratory  (MRL) 
uses  three  methods  to  carry  out  foe  analysis;  these  are  Bruceton,  Probit  and  Army  Materiel 
Command  Regulation  (AMCR). 

Much  of  the  evaluation  of  high  explosive  sensitivity,  fuse  tests,  etc  is  based  on 
collecting  a  set  of  data  that  is  obtained  at  various  levels  for  a  given  test  condition.  The  results 
of  the  individual  tests  are  recorded  as  'detonation,  non-detonation'  or  'fire,  no-fire'.  Reducing 
the  output  from  the  tests  to  a  simple  'go,  no-go'  criterion  allows  the  data  to  be  analysed 
statistically  by  the  models  cited  above.  These  three  models  are  widely  used  at  MRL  for  foe 
assessment  of  material  and  component  sensitivity.  Manually  performing  these  calculations, 
is  time  consuming  and  prone  to  human  errors,  thus  producing  misleading  results.  For  these 
reasons,  a  user-friendly  computer  model  was  developed  and  validated  against  several 
sensitivity  tests.  The  report  describes  the  software  and  it's  application  to  four  types  of 
explosive  sensitivity  tests. 


2.  Description  of  Equations  Used  in  the  Statistical 

Analysis 


The  equations  described  below  are  the  basis  of  foe  computer  model  used  to  undertake 
foe  statistical  evaluation  of  a  given  set  of  data.  All  equations  and  conditional-trees  have  been 
programmed  and  reference  graphs  digitised,  so  no  manual  calculations  are  required. 

2.1  Bruceton  Method 

The  basis  of  foe  Bruceton  Analysis  method  used  in  the  program  is  extracted  from 
Dixon  and  Mood  [  1  ].  Assuming  that  foe  data  under  analysis  is  normally  distributed,  this 
method  provides  a  good  estimate  of  foe  mean  (m)  ,  at  foe  mean  there  is  an  equal  chance  that 
the  charge  will  detonate  or  non -detonate  The  spacing  (d)  between  the  test  levels  in  the  Up- 
Down  scheme  must  be  constant  and,  ideally,  the  analysis  requires  a  large  sample  size,  say  at 
least  fifty,  to  avoid  misleading  results.  However,  in  practice,  due  to  experimental  time  and 


cost  constraints  testing  may  involve  as  few  as  ten  shots  although  twenty  to  twenty-five  shots 
is  common.  The  following  expressions  are  the  basis  of  the  program's  computations. 


A.  Calculation  of  the  mean  estimate  (m). 


m 


c  +  d 


(1) 


where, 

c,  is  the  lowest  level  used  in  Bruceton  Analysis 

d,  is  the  constant  step  interval  between  levels 

i,  is  a  number  given  to  each  test  level  starting  with  zero 
^  n, ,  represents  the  smallest  total  number  of  detonations  or  non-detonations 
±  depending  on  the  stepping  method,  a  plus  or  minus  sign  is  used. 


A  decision  to  use  either  the  plus  or  the  minus  sign  in  equation  (1)  is  based  on  the  following; 


In  a  "Increment  when  Detonate  "  stepping  method,  if  the  detonations  are  the  smaller 
of  the  two,  apply  the  plus  sign  to  the  equation  or  conversely,  the  minus  sign  is  used  when 
analysis  is  based  on  non-detonations 


In  a  "Decrement  when  Detonate  "  situation,  apply  the  minus  sign  to  the  equation 
when  analysis  is  based  on  detonations,  otherwise,  use  the  plus  sign  when  non-detonations  are 
taken  for  the  calculation 


The  mean  and  standard  deviation  are  estimated  either  from  the  detonations  or  the 
non-detonations,  whichever  is  die  smaller. 

B.  Calculation  of  the  standard  deviation  estimate  (a). 

(i)  compute  the  value  of  M, 


(ii)  M  is  then  used  to  find  the  value  of  s  from  Appendix  II  Graph  I  and  II, 
if  M  >  0.3  then,  use  the  following  identity 

s  =  1 .62  (M  +  0.029),  which  is  an  estimate  for  Graph  I  (3) 

ifA/£  0.3  then 

s  is  read  directly  from  Graph  IL,  which  has  been  digitised  and  incorporated  into  the 
program. 

If  M  falls  below  0.3,  the  curve  that  is  closest  to  diff  (eq.  3a)  is  chosen  and  the  value 
of  s  is  read  straight  from  the  curve  with  respect  to  the  M  value.  The  interpolation  method 
(eq.  3b)  is  applied  if  the  calculated  M  value  lies  between  two  digitised  M  values,  giving  a 
more  accurate  value  of  s. 


In  deciding  which  curve  on  Graph  11  to  use  to  obtain  a  value  for  s,  the  following 
procedure  is  performed.  This  is  based  on  the  equation; 

diff  =  Cnettrat  ~  m  (3a) 

d 

where,  cnearest  is  the  nearest  test  level  to  m  and  diff  represents  the  distance  between  the  mean 
and  the  nearest  level  to  the  mean.  If  diff  lies  on  one  of  the  four  curves,  the  s  value  is  read 
directly  from  that  curve,  otherwise,  the  curve  that  is  nearest  to  diff  is  chosen  to  obtain  the 
value  of  r  with  respect  to  the  value  of  M.  This  procedure  is  depicted  below,  where  h  is  one  of 
the  testing  levels. 

if  diff't  0  and  diff<  0.05 
select  m  -  h  =  0  curve, 
else 

if  diff  >  0.05  and  diff  <  0.175 
select  m  -  h  =  O  ld  curve, 
else 

if  diff>  0. 175  and  diffZ  0.375 
select  m  -  h  *=  025d  curve, 
else 

select  m  -  h  =  O.Sd  curve 

Linear  interpolation  (eq.  3b)  has  been  implemented  in  the  software  to  interpolate 
between  points  on  the  closely  digitised  curves. 

Example, 

(s2 -Si)  (3b) 

n2  ~ni 

where, 

s  is  the  interpolated  value  corresponding  to  n 
n  is  the  point  that  lie  between  Rj  and  on  an  axis 

jj,  s2  is  the  two  points  that  lie  on  the  other  axis. 

(iii)  The  standard  deviation  is  computed  as, 

a  =  ds  (4) 

C.  Calculation  of  the  mean  standard  error  (aro). 

The  mean  (m)  presented  represents  only  the  mean  of  the  sample  tested  and  not  the 
mean  of  the  whole  population,  subsequently,  it  is  exposed  to  sampling  error.  To  adjust  for 
this,  a  standard  error  for  m  is  obtained  using 


where  G  factor  is  obtained  from  Appendix  II  Graph  HI  and  IV,  with  respect  to  the  s  value. 
Depending  on  the  s  value,  the  G  factor  is  either  read  directly  from  the  digitised  points  or  from 
interpolation  between  two  digitised  G  values.  For  example,  if  an  estimated  mean  value  lies 
on  or  half  of  d  away  from  the  nearest  testing  level,  the  G  value  is  read  directly  from  the 
digitised  curve,  otherwise,  the  G  value  is  obtained  by  interpolating  between  the  curves. 
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D.  Calculation  of  95%  confidence  limit  of  mean. 

The  95%  confidence  limit  represents  a  0.95  probability  where  the  true  mean  lies  and 
it  is  a  way  to  define  the  accuracy  of  the  mean  acquired. 


L95%(m)  =  m±t095 


£"1+1.2 

I". 


where  the  value  of  f0  95  is  given  by  the  t  distribution  for  rtj  -1  degrees  of  freedom  and 

r  y«i +1.2] 

== -  is  a  correction  factor  which  is  ineffective  when  large  sample  size  are  used. 


2.2  Probit  Method 

The  basis  of  die  equations  implemented  in  die  software  is  originated  from  Finney  [2]. 
Transformation  of  percentages  to  probit  is  the  basis  of  this  statistical  analysis  method.  It 
assumes  that  data  samples  are  either  normal  or  log-normal  distribution.  The  method  is 
capable  of  being  adapted  to  data  of  unequal  sizes  at  each  test  level  and  the  difference  between 
these  levels  need  not  be  equally  spaced 

For  a  given  percentage,  the  probit  (  Y)  may  be  defined  as 


where. 


x  is  the  test  level,  ie.  voltage  level  for  cap  sensitivity  test 
pp  is  the  mean  of  the  population 

Op  is  the  standard  deviation  of  the  population 


(7) 


I  *1*0 

The  second  part  of  equation  (7),  - —  ,  represents  the  percentage  of  the  area  under  the 

L  aP 

normal  distribution  curve  for  various  x  values  and  is  commonly  written  as  ap  +ppx,  where 


subscript  'p'  denotes  the  population. 
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A.  Calculation  of  the  Empirical  Probit. 

The  transformation  of  observed  percentages  to  probit  is  designated  as  Empirical 
Probit.  Using  Abramowitz  and  Stegun  (3],  which  approximates  the  inverse  of  the  normal 
distribution  function,  the  following  applies; 


i 

*  JL 


h  \‘  ldy 


where,  y  = 


Approximation  of  the  function  from  Abramowitz  and  Stegun  [3,  eq.  26.2.23]  is 
derived  as, 

aLtaf.^.J  (<>) 

l  +  dgf+djf  +^1 

where. 


1=1-  2  In 


1 2  2) 


%  observed  percentage  at  each  test  level 

c0  2.515517 

c,  0.802853 

c2  0.010328 

d0  1.432788 

dx  0.189269 

d2  0.001308 

Therefore,  T  =  5.0  +  -ep,  when  0  <  %  <  0.5 
K  =  5.0  +  ep,  when  0.5  <%<  1.0 

B.  Calculation  of  the  Provisional  Probit. 

The  provisional  line  provides  a  means  of  estimating  the  level  at  which  a  specified 

Expected  Probit  occurs.  The  equation  for  the  provisional  line  is  computed  using. 


where, 


T=5.0+y 

y  =  V(x-vJ: 

P  =  —  ;  is  the  slope  of  the  provisional  line 
a 

p  is  the  sample  mean  estimate 
o  is  the  sample  standard  deviation  estimate 
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C.  Calculation  of  the  Expected  Probit. 

The  Expected  Probit  for  each  empirical  percentage  versus  its  Provisional  Probit  is 
computed  as 

r  =  5.0+/  (II) 

where,  y"  =/  +  %  -- 
zy 

%  observed  percentage 

p  is  defined  in  equation  (17) 

Zy  is  defined  in  equation  (19) 

D.  Calculation  of  the  threshold  statistics  for  a  given  percentage  (x*). 

The  mean  of  the  no-fire  level  is  calculated  using  equation  (12),  and  is  originated  from 
equation  (8a).  The  value  for  ep  is  derived  with  equation  (9),  which  is  equivalent  to  the  no¬ 
fire  probability. 

x%  =  n  +  epa  (12) 


E.  Calculation  for  the  threshold  standard  error  for  a  given  percentage  (  am). 


This  is  given  by. 


l  I  1  |  (*% 

MX*/  s» 


(13) 


where  w/is  the  weight  factor  as  in  equation  ( 16);  ntwc 

x  is  defined  in  equation  (20);  weighted  mean  of  x 
Sxx  is  expressed  in  equation  (22) 

F.  Calculation  of  the  95%  confidence  limit  of  the  threshold  (  x95%). 

The  95%  confidence  limit  of  the  no-fire  percentage  threshold  is  calculated  using. 


x95  %=x%-t09Sam  (14) 

Assuming  large  sample  theory,  the  value  of  r0.9S  front  the  t  distribution  table  is  equal 
to  1.645.  For  the  95  percent  confidence  limit  of  the  fire  threshold,  replace  equation  14  with  a 
plus  sign. 

2.2.1  Computational  Procedures. 

1.  The  Empirical  Probit  for  each  test  level  with  respect  to  the  no-fire  percentage  is  calculated 
using  equation  (9). 

2.  The  Provisional  Probit,  the  Expected  Probit  together  with  the  no-fire  threshold,  the 
variance  and  the  95%  confidence  limit  are  calculated  as  follows, 

2.1  get  current  estimates  (it,o), 

2.2  calculate  Provisional  Probit  using  equation  ( 10), 

2.3  evaluate  for  the  weight  coefficient  (  wc)  and  weight  factor  (wy)  for  each 
Provisional  Probit  value  as  follows. 


(15) 
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ntwc 


(16) 


wf 


where,  q  =  1  -  p,  and  n,  is  the  total  number  of  shot  at  each  level 

The  value  of  p  in  equation  (IS)  is  derived  as, 

ify<0,  then  p-f 
ify  £  0,  then  p  =  \  -f 


where,  /is  the  Normal  Probability  Density  Function,  and  is  defined  as 


(17) 


f  =  Zy  [uo/r+ax/r2  +«2A3]  in  A&S[3,  eq.  26.2.16] 


where. 


a  0.33267 

a0  0.4361836 

ax  -0.1201676 

a2  0.937298 


l+«|y| 


(18) 


(19) 


2.4  the  Expected  Probit  for  each  empirical  percentage  versus  its  Provisional  Probit  is 
evaluated  using  equation  (1 1). 

2.5  compute  for  probit  line  and  the  sum  of  squares  and  products  about  the  mean. 


'y.wfx 

x  =  ^  ;  weighted  mean  of  x 

ZuWf 

(20) 

y  =  *=  —  ; weighted  mean  of  y 

2>/ 

(21) 

a.  =  xy  -  {^zf 
^  ’  I-/ 

(22) 

%  - 

2>/ 

(23) 

c  ..2  (Xw/y")z 

Syy  ~  Zuw(y  ~ 

2>/ 

(24) 
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2.6  calculation  for  the  new  estimate  (  p,o)  and  chi-square. 


i  - 

P  Sxy 


H  =  x  -oy 

..  (Sxy)2 

chi  -  square  =  Syy - 

Sxx 


(25) 

(26) 
(27) 


3.  Repeat  procedure  2. 1  to  2.6  until  the  probit  line  is  an  approximation  that  results  in 
convergence.  In  the  program,  convergence  is  reached  when  current  chi-square  minus 
previous  chi-square  is  less  than  0.0003. 

4.  Finally,  the  desired  no-fire  probability  percentage  (%  )  is  subjected  to  equation  (12),  (13) 
and  (14)  to  obtain  the  necessary  no-fire  threshold  (  x*),  the  standard  deviation  (ow)  and  the 
95%  confidence  limit  ( x95%)  for  the  threshold  statistics. 


2.3  Army  Materiel  Command  Regulation  (AMCR)  Method 

Army  Materiel  Command  Regulation  (AMCR)  [4]  is  a  simple,  fast  method  and 
capable  of  handling  large  data  at  each  test  level.  Like  Bruceton,  the  step  interval  between  test 
levels  must  be  equally  spaced.  Within  these  test  levels,  it  must  have  one  level  having  all  its 
results  going  one  way  and  another  going  the  other  way,  ie.  100%  no-fires  and  100%  fires. 

The  equations  below  are  die  basis  of  die  method  used  in  the  software. 

A.  Calculation  of  the  mean  estimate  (m). 

This  is  given  by, 

m  *  ]£  l><f  +  |^H100%  +  (28) 

where, 

y,P,  is  the  sum  of  the  no-fire  percentage  at  each  level 

d  represents  the  step  interval  between  levels 
//100%  is  the  first  test  level  that  all  sample  misfire 

B.  Calculation  of  the  Standard  Deviation  (  a). 

This  is  given  by, 

a  -  o'  m 

where,  _ 

o'  -  JlM,  -  Xp,1  ; 

is  the  standard  deviation  unoorrected  for  the  test  interval 

y  kfPj  is  the  sum  of  the  variance  factor  kt  multiplied  by  the  no-fire  percentage  (PJ  at  that 
level.  At  the  all  no-fires  test  level,  kj=l  and  k^kj+2  for  subsequent  levels,  ie.  it15 1, 3,5,7,.... 
y  Pj  is  the  sum  of  the  no-fire  percentage  at  die  current  level. 


(29) 

(30) 
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2.4  Choice  of  Method 

The  three  methods  described  above  provide  a  means  of  estimating  and  evaluating 
data.  However,  one  may  be  better  suited  for  a  particular  set  of  data  than  another.  Here  are  a 
few  points  which  might  be  helpful  in  deciding  which  method  to  apply; 

(i)  Bruceton  Analysis  : 

•  concentrates  very  much  on  the  mean,  hence  provides  more  accurate  estimate  about 
the  mean.  But  it  is  not  a  good  method  for  estimating  small  or  large  percentage  points 

M- 

•  at  least  fifty  data  points[  I  ] . 

•  requires  fewer  tests  for  a  given  accuracy!  1  ]. 

•  the  95%  confidence  limit  is  implemented. 

•  necessary  to  have  constant  spacing  between  test  levels;  software  restriction. 

•  maximum  number  of  data  points  is  two  hundred,  but  expandable;  software 
restriction. 

(ii)  Probit  Analysis: 

•  does  not  require  constant  spacing  between  test  levels  [2]. 

•  many  data  points  at  each  test  level  is  possible. 

•  maximum  thirty  test  levels,  but  expandable;  software  restriction. 

•  the  95%  confidence  limit  is  implemented. 

•  can  not  utilise  the  all  fires  and  all  no-fires  levels.  Fisher  and  Yates  (1964,  Table  IX). 

(iii)  AMCR  Analysis: 

•  constant  spacing  between  test  levels  is  necessary  [4], 

•  thirty  test  levels  maximum,  but  expandable;  software  restriction. 

•  quick  and  easy  way  to  approximate  results. 

•  provides  estimate  of  the  mean  and  standard  deviation  only  [4], 

•  has  no  95%  confidence  limit. 

•  requires  the  all  fires  and  all  no-fires  levels  [4]. 

3.  Description  of  the  Statistical  Analysis  Software 

3.1  Software  and  Hardware  Requirements 

The  recommended  computer  system  should  include  : 

-  IBM  AT  or  compatible  Personal  Computer 

-  IBM  colour  VGA  graphics  card  or  any  other  card  compatible  with  Microsoft 
Windows,  version  3.0  or  higher 

-  IBM  colour  monitor  or  compatible.  Under  Windows  environment  it  should  be 

running  at  VGA  mode  (640  x  480  pixel) 

-  at  least  2  megabytes  of  Random  Access  Memory 

-  Hard  disk  with  at  least  1  megabyte  available 

-  Microsoft  mouse  or  any  other  pointing  device  compatible  with  Microsoft 
Windows  3.0  or  higher 

-  HP  LaserJet  II  printer  or  higher  or  any  other  compatible  LaserJet  printer 
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With  Operating  System  (PC  or  MS-DOS  3.0)  and  Microsoft  Windows  version  3.0  or 

3.1  installed  on  the  computer  system,  the  software  package  requires  two  more  files;  the 
executable  file  (stwin93.exe)  and  Borland's  Dynamically  Linked  Library  file  (bwcc.dll). 

3.2  General  Development  Details 

The  three  methods  of  statistical  data  analysis  (Bruceton,  Probit  &  AMCR)  were 
converted  into  a  user-friendly  computer  program.  Two  versions  of  the  program  were 
developed;  a  DOS  version  and  a  Windows  version.  However,  the  following  descriptions 
apply  to  the  Windows  version  only.  The  Windows  version  of  the  program  is  very  easy  to 
use.  It  has  been  written  to  the  Windows  standards,  meaning  that  die  users  deal  with  the  same 
screens,  menus  and  methods  regardless  of  the  application  they're  working  in,  thus  reducing 
the  teaming  curve  and  boosting  productivity. 

The  application  program  was  developed  under  Microsoft  Windows  version  3.0 
environment  using  Borland  C++  version  3. 1  software  development  tools.  During  the 
development  phase,  the  program  was  compiled  and  run  on  a  486  computer  system  with  a 
super  VGA  monitor  running  at  640  x  480  pixel  as  well  as  eight  megabytes  of  Random  Access 
Memory.  The  Output  device  used  was  a  HP  LaserJet  III  printer.  Of  course,  this  is  not  the 
minimum  criterion.  Refer  to  section  3. 1  for  software  and  hardware  requirements. 

The  computer  model  requires  two  files  to  run,  the  executable  file  (stwin93.exe)  and 
Borland's  Dynamically  Linked  library  file  (bwcc.dll).  During  the  development  of  die 
application  software,  interaction  with  operators  enabled  the  refinement  of  the  software,  for 
example,  having  the  up-down  graph  or  ratio  bar  graph  helper,  nake  the  raw  data  more 
readable.  The  program  has  been  tested  and  is  now  being  used  in  various  analyses  at  MRL. 
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3.3  Application  Software  Menu  Descriptions 

When  the  statistical  analysis  software  is  loaded,  the  main  application  window  will 
appear  on  the  screen,  as  shown  in  Figure  1 .  The  program's  main  menu  bar  has  four  menu 
items  (File,  Method,  Misc. .  and  Help)  that  are  available  to  the  operator.  Each  menu  item 
contains  more  sub-menu  items  to  perform  a  specific  task.  This  is  depicted  in  Figure  2  and  an 
explanation  of  each  menu  command  is  described  below.  The  structured  diagrams  of  the 
program  can  be  referred  to  in  Appendix  I. 


MflL  Statistical  Data  Analysis 


file  Method  Miss...  Help 


■s 


Material  Research  Laboratory 

Statistical  Analysis  Software 


Figure  1.  Main  application  window  with  four  menu  items 
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Eruceton  Title  Author 
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AMCR  Custom 

Terminology 
Convert 
Convert  getUp 
Standard 
Custom 


Figure  2.  Program's  menu  items 
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3 3.1  File  Menu 


File  menu  enables  users  to  toad  existing  data  files,  save  data,  select  an  output  device  for 
printing,  and  exit  the  application  program. 

File)  Open  command  enables  the  user  to  quickly  recall  a  previously  saved  data  file  from 
disk.  If  the  required  data  file  is  not  listed  on  the  Files  list  box,  select  a  new  drive, 
directory  or  select  a  different  file  type. 


Figure  3.  File  Open  dialog  box 


Filename  input  box  prompts  the  user  to  enter  die  filename  explicitly  or  using  the  DOS 
wildcards  such  as  *  and  ?  to  filter  out  unwanted  files.  For  more  information  on  DOS 
wildcards,  refer  to  DOS  documentation. 

Flies  list  box  lists  files  in  die  current  directory.  Depending  on  what  is  on  die  input  box, 
various  filenames  will  appear  in  this  list  box.  When  die  program  is  initially  called  up,  the 
wildcard  appear  on  the  input  box  is  *.*  which  will  list  all  filenames  in  current  directory. 
Once  an  analysis  method  is  selected,  die  filename  extension  is  changed  to  an  appropriate 
one,  ie.  *.BRU  for  Bruceton,  ’.PRO  for  Probit  and  *.AMC  for  AMCR  file  type.  Hence, 
only  die  relevant  file  type  is  displayed. 

Directories  list  box  lists  the  directories  and  drives  available  on  the  system. 

Note:  In  opening  an  existing  data  file,  the  default  or  current  settings  will  be  changed  to 
the  settings  specified  in  that  file.  This  is  extremely  handy  in  terms  of  editing  the 
Experimental  dialog  box  (Figure  9). 

File) Save  command  saves  the  current  settings,  experimental  information  and  raw  data  to 
disk  under  the  specified  filename.  The  first  time  the  data  is  saved,  the  Save  command 
function  as  the  SaveAs  command.  Refer  to  Appendix  III  for  file  format  details. 

For  easy  file  referencing,  it  is  recommended  the  following  file  extensions  should  be  used. 


18 


Bruceton  data  file,  filename.BRU 

Probit  data  file,  filename.PRO 

AMCR  data  file,  filename.AMC 

File] Save  As  command  enables  the  user  to  save  an  unnamed  data  file,  ie.  new  data  set,  or 
save  an  existing  or  updated  data  files  under  a  different  name,  different  directory  or 
different  disk  drive.  In  the  input  box,  type  in  the  filename,  select  die  disk  drive  and 
directory  required,  then  select  OK. 


Figure  4.  SaveAs  dialog  box 


Filename  input  box  prompts  the  user  to  enter  a  new  filename. 

Directories  list  box  lists  die  directories  and  drives  available  on  the  system  in  which  data 
may  be  saved. 

Filej  Print  command  enables  the  user  to  produce  hard  copies  of  analysed  results.  To  cancel 
current  print  job,  select  Cancel  button  in  the  Printing  dialog  box  that  will  appear  on  the 
screen  until  the  current  print  job  is  finished  printing.  It  is  a  good  practise  to  save  the  data 
before  printing.  That  way,  if  a  printer  error  or  other  problems  occur,  data  is  not  lost. 
Copies  of  the  printed  output  from  the  HP  LaserJet  III  printer  are  shown  in  Appendix  IV. 

FilejPrinterSet  Up  command  enables  the  user  to  select  a  printer  that  has  already  been 
installed  on  die  system.  For  more  information  on  installing  printer  driver  software,  refer 
to  Windows  documentation.  When  selected,  it  calls  up  a  Select-Printer  dialog  box 
showing  the  current  default  printer  used  as  in  Figure  5 . 


Figure  5.  Select-Printer  dialog  box 


Down  Arrow  (^)  button  lists  the  available  printers)  available  for  use.  Here  foe  user  can 
select  a  preferred  printer  as  a  default  printer. 

Setup  button  calls  up  a  separate  Printer  Setup  dialog  box  which  enables  foe  user  to 
modify  printer  settings,  such  as  different  paper  size,  orientation,  source  of  paper,  graphics 
resolution,  etc.  For  more  information,  refer  to  Windows  documentation.  The  settings 
used  to  obtain  foe  printed  copies  in  Appendix  IV  are  as  follows. 
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Figure  5a.  Printer  Setup  dialog  box 


FikjExit  command  quits  the  application  program  and  returns  to  the  previous  active 
window.  When  the  exit  command  is  selected,  a  Confirmation  dialog  box  appears  on  the 
screen,  to  confirm  the  exit  command.  Furthermore,  if  the  data  is  not  saved,  it  will  prompt 
the  user  to  save  the  data. 


3.3.2  Method  Menu 

Method  menu  enables  the  user  to  select  the  three  methods  of  analysis  available. 

Method]  Bruceton  command  enables  the  user  to  invoke  the  Bruceton  analysis  method. 
When  selected,  a  Bruceton  input  dialog  box  appears,  allowing  the  user  to  enter  analysis 
conditions  and  the  data  to  be  processed.  (Figure  6) 


Figure  6.  Bruceton  input  dialog  box 

Increment  or  Decrement  radio  button  indicates  the  type  of  stepping  method  used  in  die 
Bruceton  up-down  scheme.  For  example.  Increment  checked  means  step  up  a  level 
when  detonation  occurs,  or  alternatively.  Decrement  checked  means  step  down  a  level 
when  detonation  occurs.  Default  is  Increment  checked. 

Step  Interval  input  box  requires  the  user  to  enter  the  stepping  size  between  levels. 

First  Test  Level  input  box  requires  the  user  to  enter  the  first  Bruceton  test  level.  The 
combination  of  these  two  inputs  permits  the  automatic  data  determination  feature  to  be 
implemented. 

Shot  Detail  input  box  is  optional.  This  feature  is  implemented  to  allow  cross  reference 
between  the  input  data  and  the  experimental  records.  Three  recognisable  formats  are 
implemented; 

First  is  the  default,  where  there  is  no  value  entered  in  the  input  box,  it  assumes  the 
first  shot  starts  at  1  and  then  increments  by  one  in  subsequent  shots,  ie.  1,  2,  3,  etc. 

The  second  format  is  where  the  user  enters  the  initial  shot  number,  and  the  number  is 
automatically  incremented  by  one  between  shots,  ie.  168,  169,  170,  etc. 
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The  third  format  permits  character  input,  where  the  user  supplies  the  initial  shot 
index.  This  shot  index  must  be  preceded  by  two  characters,  ie.  SC10001,  SC10002, 

SC  10003,  SC  1004,  etc 

Of  course,  the  three  formats  can  be  mixed  together,  ie.  168,  169,  2,  3,  SC1001,  170, 
SC1002,  etc.  All  three  formats  have  a  maximum  of  eight  characters. 

Fire  button  implies  that  a  shot  resulted  in  detonation,  and  the  level  increments  or 
decrements  accordingly. 

NoFire  button  implies  that  a  shot  resulted  in  non-detonation,  and  again,  the  level 
increments  or  decrements  accordingly. 

Bruceton  level  list  box  lists  the  shot  number,  shot  details,  test  level  and  result  of  the 
shots  recorded.  A  scroll  bar  will  appear  on  the  right  side  if  the  number  of  shots  recorded 
exceeds  die  length  of  the  list  box.  11181  way,  the  user  can  scroll  up  or  down  the  long  list 
of  shots. 

To  alter  any  recorded  shot,  the  user  double-clicks  the  mouse  on  the  required  shot. 
This  produces  a  separate  dialog  box  asking  for  Delete  or  Modify  command.  The  Modify 
command  modifies  the  result  of  die  shot  and  automatically  updates  any  shots  that 
follow.  The  Delete  command  can  be  earned  out  only  if  die  but  recorded  shot  is 
selected. 

Note:  A  limit  of  two  hundred  shots  has  been  implemented. 

OK  button  closes  die  Bruceton  input  dialog  box,  saves  the  data  to  memory,  performs  the 
calculation,  and  displays  the  results  and  input  parameters  onto  a  Bruceton  result 
window.  A  sample  of  Bruceton  result  window  is  shown  in  Figure  6a. 

Caned  button  closes  the  Bruceton  input  dialog  box  and  returns  control  to  die  main 
application  window. 

Note:  No  inputs  or  settings  can  be  altered  in  the  middle  of  data  input,  ie.  die  step 
interval,  die  first  test  level  input  box  and  die  Increment/Decrement  radio  button  cannot 
be  altered  after  die  Fire  or  NoFire  button  is  depressed. 
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Bruceton  Results-  '7^  1 
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Shot  #  S -detail  .Level!: Result Ifc* 
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Figure  6a.  Bruceton  result  window 


Add  button  permits  the  user  to  add  ranging  shots  to  existing  Bruceton  data. 

Remove  button  enables  the  user  to  delete  existing  ranging  shots. 

Command|  GoBack  command  in  this  window  allows  the  user  to  go  back  to  Bruceton 
input  dialog  box  to  correct  or  update  data,  and  returns  with  an  updated  Bruceton  result 
window 

Command|  Done  quits  the  Bruceton  result  window  and  returns  to  the  main  application 
window. 

Note:  The  ranging  shot(s)  is  for  user  information  only,  it  is  not  processed  by  the 
Bruceton  analysis.  On  the  hardcopy,  the  ranging  shots  are  distinguished  from  the 
Bruceton  data  by  an  asterisk  (*)  preceding  the  shot  number. 
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Method)  Probit  command  enables  the  user  to  invoke  the  probit  analysis  method.  When 
selected,  a  probit  input  dialog  box  appears,  allowing  the  user  to  input  probit  data. 


Figure  7.  Probit  input  dialog  box 

Log  or  NonLog  radio  button  enables  die  user  to  select  either  log  or  non-log  test  levels 
before  the  probit  calculation  begins.  To  avoid  negative  results  a  log-normal  distribution 
is  recommended.  Therefore,  Log  is  set  as  default. 

Note:  If  the  samples  are  already  log-transformed,  select  the  NonLog  setting. 

Add  button  lets  die  user  add  data  to  die  probit  list  box.  A  maximum  of  thirty  levels  is 
permitted.  Depressing  the  Add  button  calls  up  another  probit  formatted  dialog  box  (see 
Figure  7a),  which  requires  the  user  to  enter  die  value  for  the  test  level,  the  total  number 
of  fires  and  no-fires  at  each  level.  OK  button  closes  the  probit  formatted  dialog  box,  and 
adds  the  data  to  die  probit  list  box.  Caned  button  closes  the  probit  formatted  dialog  box, 
without  adding  data  to  the  probit  list  box. 


Figure  7a.  Probit  formatted  dialog  box 


Delete  button  lets  the  user  remove  recorded  data  from  the  probit  list  box.  To  remove 
data  from  die  list,  select  the  data  to  be  deleted  using  the  mouse  and  then  depress  the 
Delete  button. 
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OK  button  closes  the  probit  input  dialog  box,  saves  and  processes  the  data  using  the 
probit  analysis  method.  Probit  result  window;  depicted  Figure  7b.,  displays  various 
calculated  parameters  and  results. 

Cancel  button  closes  the  probit  input  dialog  box  and  returns  to  the  main  application 
window. 


Figure  7b.  Probit  result  window 

Command|  GoBack  and  Command)  Done  functions  exactly  the  same  as  described  in 
Bruceton  result  window  section. 

Note:  To  modify  data,  double-click  the  data  point  using  the  mouse.  The  probit 
formatted  dialog  box  will  appear  allowing  the  user  to  perform  the  necessary  changes. 
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Method]  AMCR  command  enables  the  user  to  invoke  the  AMCR  analysis  method. 
When  selected,  an  AMCR  input  dialog  box  appears,  allowing  the  user  to  enter  AMCR 
data  for  processing 


The  AMCR  input  dialog  box  has  a  similar  layout  and  functions  as  that  of  the  probit 
input  dialog  box,  except  that  the  Log  or  NonLog  radio  button  does  not  applied.  For 
more  information,  refer  to  probit  command  section.  The  AMCR  result  window  is  shown 
in  Figure  8.  The  necessity  to  have  the  first  test  level  all  no-fires  and  the  last  test  level 
all  fires  is  depicted  below. 


Figure  8.  AMCR  result  window 


33.3  Misc...  Menu 

Miscellaneous  menu  enables  the  user  to  set  the  preferred  terminology,  input  experimental 
details,  and  update  old  file  format. 

Misc..  .|Title  lets  the  user  enter  experimental  details  to  be  included  as  part  of  the  printout. 
There  are  two  experimental  detail  formats;  Standard  and  Custom  format. 

Misc... (Titlej Standard  command  calls  up  a  standard  experimental  detail  dialog  box.  This 
dialog  box;  as  in  Figure  9,  has  the  standard  field  names  such  as  explosive  composition, 
density,  etc.  All  the  user  needs  to  do  is  to  type  in  the  relevant  information 
corresponding  to  field  names.  Of  course,  this  is  optional.  If  no  information  is  entered, 
no  detail  is  printed  after  the  field  names. 
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Figure  9.  Standard  experimental  detail  dialog  box 

Clear  button  enables  the  user  to  quickly  clear  all  information  appearing  on  the  input 
boxes.  Once  cleared,  all  information  is  erased  from  the  memory.  Therefore,  for  minor 
alterations  edit  the  field  information  individually. 

Note:  Experimental  details  are  stored  with  each  data  file  and  can  be  recalled. 

OK  button  closes  the  dialog  box  and  saves  all  the  information  that's  on  the  input  boxes. 

Cancel  button  closes  the  dialog  box  without  saving  the  information. 

Misc...|Title| Custom  command  calls  up  a  custom  experimental  detail  dialog  box.  In  this 
empty  dialog  box,  the  user  needs  to  enter  the  preferred  field  names  as  well  as  relevant 
information  with  respect  to  f  -  field  names.  Refer  to  Appendix  IV,  Figures  4,5,6, 7  and  8 
for  examples.  Again,  if  no  field  name  or  information  is  entered,  a  blank  space  is  printed. 
Buttons  that  appear  on  this  dialog  box  has  the  same  function  as  that  of 
Misc . .  |T itlel  Standard . 


Misc... | Preference  command  enables  the  user  to  select  predefined  terminology  and  the 
units  system  to  appear  on  the  printout.  When  selected,  the  dialog  box  allows  the  user  to 
choose  the  preferred  system. 


Figure  10.  Preference  dialog  box 

Predefined  Terminologies: 

Detonation  and  NonDetonation  -  D  and  ND  *  | 

Fire  and  NoFire 

FandNF 

Hit  and  Miss 

H  and  M 

Success  and  Failure 

S  and  F 

Predefined  Unit  system  : 

None 

Blank  * 

Millimetre 

mm 

Centimetre 

cm 

Meter 

m 

Inch 

in 

Decimal  Precision : 

Thousandth 

3  decimal  places  * 

Hundredth 

2  decimal  places 

Tenth 

1  decimal  places 

Zero 

No  decimal  place 

*  denotes  default  settings 

1 

Misc...|Convert  command  enables  the  user  to  convert  an  existing  old  data  file  format  to 
current  data  file  format.  There  are  two  forms  of  data  file  conversion;  Standard  and 

Custom  format.  The  difference  between  the  two  is  that  Standard  conversion  assumes  all 
settings  use  the  default  value  (standard  field  names,  Detonation/NonDetonation,  no  unit 
system,  etc.),  whereas.  Custom  conversion  requires  the  user  to  manually  fix  the  settings 
through  various  dialog  boxes  during  the  conversion  process.  These  dialog  boxes  are  the 
preference  select  and  the  experimental  detail  input  dialog  boxes  Refer  to  above  for 

dialog  box  details. 
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Note:  All  data  files  saved  using  the  DOS  version  of  the  computer  model  are  considered 
old  file  format. 

Misc...|Convert  SetUp  command  enables  the  user  to  select  the  conversion  process. 

These  are  the  type  of  data  file  and  the  type  of  format  conversion.  When  selected  it  calls 
up  a  Con vertS ctup  dialog  box.  Here  the  user  can  select  die  type  of  data  file  (Bruceton, 
Probit  or  AMCR  type)  and  the  format  of  the  converted  file  (Standard  or  Custom). 
Default  setting  uses  Bruceton  file  type  and  the  Standard  format. 

33.4  Help  Menu 

Help  menu  enables  the  user  to  call  up  help  information. 

HelpjAuthor  command  calls  up  a  programmer  detail  dialog  box  for  enquiries. 

3.4  Limitations  and  Extensions 

In  order  to  obtain  reliable  Bruceton  results,  at  least  fifty  data  points  are  ideally 
required.  The  computer  model  permits  up  to  two  hundred  data  points,  which  should  be 
sufficient  for  most  tests.  Likewise,  a  maximum  of  thirty  ranging  shots  are  possible.  Due  to 
space  limitations  on  A4  size  paper,  a  maximum  of  fifty  shots  can  be  printed  in  the  table,  but 
the  whole  Bruceton  series  is  graphically  represented. 

Probit  and  AMCR  data  analysis  have  a  limit  of  thirty  levels,  and  die  total  number  of 
data  points  on  each  level  should '  less  than  one  thousand.  A  minimum  of  three  levels  are 
required. 

Note:  All  limitations  above  are  *  .itware  restriction  and  are  extendable. 

Further  extensions  to  the  computer  model  could  incorporate  other  commonly  used 
statistical  data  analysis  methods  and  the  enhancement  of  existing  algorithms. 


4.  Application  of  the  Computer  Model  to  Explosive 

Testing 


The  software  package  has  been  successfully  applied  to  several  tests  at  MRL 
concerned  with  estimating  the  sensitivity  of  explosive  materials  and  the  reliability  of 
explosive  devices.  As  a  consequence  of  the  user  evaluation,  the  software  was  developed  to 
Windows  standards,  so  it  is  easy  to  use,  yet  flexible  enough  to  accomondate  various  tests. 
And  the  printout  format  has  been  tailored  to  produce  hardcopy  in  a  form  ready  for  reporting, 
presentation  and  storage.  The  tests  serviced  by  the  statistical  software  package  are  described 
below  and  examples  of  the  resultant  printout  are  given  in  Appendix  IV. 

4.1  The  Jet  Sensitivity  test 

4.1.1  General  description  of  test. 

The  jet  sensitivity  test  [S]  was  developed  at  MRL  to  investigate  the  sensitivity  of 
both  bare  and  covered  explosives  to  shape  charge  jets.  The  experimental  set-ups  for 
assessing  the  sensitivity  of  both  covered  and  bare  explosives  are  depicted  in  Appendix  V, 
Figure  1  and  2  respectively. 
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The  covered  explosive  configuration  was  designed  to  simulate  the  cross  section  of  a 
munition  case  and  explosive  filling  penetrated  by  a  shape  charge  jet.  The  confinement  of  the 
near  and  far  side  of  the  munition  casing  is  represented  by  the  cover  and  witness  plate 
respectively.  The  standoff  distance  between  the  shape  charge  and  the  cover  is  to  allow  the  jet 
to  develop  after  the  shape  charge  is  initiated.  It  is  important  that  the  geometry  of  the  receptor 
is  of  sufficient  size  to  avoid  complications  from  shocks  reflected  from  the  side  and  witness 
block  interfaces. 

A  similar  set-up  is  applied  to  bare  explosive  configurations  except  that  an  additional 
15  mm  air  gap  is  incorporated  between  the  cover  and  the  receptor  explosive.  This 
arrangement  was  designed  to  remove  the  effect  of  the  precursor  waves  (resulting  from  jet 
impact  and  jet  penetration  of  the  barrier)  [6]  from  affecting  the  explosive. 

The  sensitivity  of  the  explosive  is  expressed  as  the  thickness  of  the  cover  material 
which  produces  a  50%  probability  of  detonation.  Variation  of  the  cover  thickness  changes 
the  jet  characteristics,  hence  the  velocity  of  the  jet  and  shock  entering  die  explosive.  Thus 
increasing  the  cover  thickness  decreases  the  jet  and  precursor  wave  velocities  and  vice  versa. 

As  in  other  tests  of  critical  'go,  no-go'  behaviour,  the  results  must  be  interpreted 
statistically.  A  decision  to  increase  or  decrease  the  cover  thickness  is  depended  on  the 
previous  test  result.  A  deep  rounded  dent  on  the  witness  plate  indicates  that  the  test 
explosive  detonated,  whereas  no  dent  indicates  a  non-detonation.  If  in  the  test  procedure,  the 
charge  detonates  then  the  cover  thickness  is  increased  by  a  fixed  thickness  and  if  it  fails  to 
detonate  then  the  cover  thickness  is  decreased  by  the  same  thickness.  Such  an  'up-down' 
procedure  is  continued  until  a  regular  pattern  of  detonations  and  non-detonations  is  obtained. 
The  results  are  then  analysed  statistically  [1]  to  give  the  mean  cover  thickness  the  standard 
error  and  the  95%  confidence  limit.  This  critical  cover  thickness  represents  the  detonation 
and  non-detonation  threshold  for  the  receptor  explosive. 

4.13.  Application  of  the  software  package. 

A  typical  test  schedule  for  a  MRL  jet  sensitivity  test  is  shown  in  Appendix  IV, 

Figures  1;  this  includes  die  results  for  the  statistical  analysis  of  the  data.  The  hardcopy  from 
the  software  of  the  results  is  shown  in  Appendix  IV,  Figures  2a  and  2b.  Figure  2a  represents 
a  set  of  data  collected  from  a  jet  sensitivity  test  on  covered  Composition-B  explosive. 

Ideally,  a  large  data  sample  is  required  to  give  a  reasonable  estimate.  However,  in  practice, 
ten  to  twenty-five  samples  are  usually  the  case.  It  uses  the  Bruceton  analysis  to  determine 
the  critical  cover  thickness  (  m),  standard  error  (Sm)  and  the  95%  confidence  limit  ( L95%(/n)). 

Appendix  IV,  Figure  2a,  illustrates  the  explosive  under  test  was  a  38  mm  diameter  by 
50  mm  length  Composition-B  cylinder  and  drat  two  cylinders  were  stacked  on  top  of  each 
other  to  give  a  total  charge  length  of  100  mm.  Hie  test  set-up  was  type  'C',  and  the  explosive 
was  covered  with  a  mild  steel  barrier. 

The  numerical  table  displays  the  detail  of  the  shots  performed  on  a  particular  test. 
Starting  with  a  63.6  mm  mild  steel  cover  thickness  (Test  Level),  the  result  was  a  non¬ 
detonation,  ND,  the  cover  thickness  was  then  decreased  by  6.4  mm  (Step  Size)  to  57.2  mm 
where  the  charge  detonated,  D.  The  subsequent  cover  thickness  was  increment  or  decrement 
accordingly  for  twenty  one  shots.  The  printout  shows  that  the  step-up  scheme  was  employed 
to  the  Bruceton  Analysis  when  a  detonation  resulted.  The  graphical  representation  provides  a 
clearer  view  of  the  data  set.  The  "Shot  Detail"  column  represents  the  user's  own  reference  or 
catalogue  number  of  each  individual  shot.  Notice,  there  are  no  ranging  shots  present  on 


figure  2a,  but  in  most  cases  ranging  shots  are  unavoidable,  as  depicted  on  Figure  2b.  There 
are  three  ranging  shots  on  Figure  2b.  They  are  shot  sequence  numbers  (No.)  16  to  18  and  are 
highlighted  by  **'. 

The  Statistical  Results  box  displays  the  mean  (m),  the  standard  error  (Sm),  and  the 
95%  confidence  limit  of  the  mean.  Other  supporting  parameters  such  as  the  standard 
deviation  (sd),  the  95%  confidence  limit,  s,  M,  and  G  values  are  also  shown.  Refer  to 
Section  8  for  symbol  definition. 

This  neatly  compacted  printout  is  customised  to  encapsulate  all  the  necessary  test 
details  and  statistical  data.  Thus,  the  printout  shows  immediately  how  the  test  was 
performed,  what  apparatus  was  required  to  carry  out  the  test  and  what  result  was  achieved. 

4.2  Gap  Test 

4.2.1  General  description  of  test. 

The  MRL  Small  Scale  Gap  Test  (SSGT)  provides  a  simple  way  to  determine  the 
sensitivity  of  an  explosive  to  initiation  by  shock  waves.  A  donor  charge  generates  a  shock 
pressure  of  uniform  magnitude  that  is  transmitted  to  an  acceptor  charge  through  an 
attenuating  barrier.  Commonly,  die  inert  material  for  the  barrier  consists  of  laminated  brass 
shim.  Variation  to  the  thickness  of  this  barrier  alters  the  strength  of  the  shock  waves  entering 
the  acceptor,  hence  the  result  of  the  test.  In  doing  so,  one  can  determine  the  critical  barrier 
thickness  required  to  initiate  detonation  in  acceptor  charge.  The  larger  the  critical  gap,  die 
more  shock  sensitivity  is  the  acceptor  material. 

To  determine  an  accurate  value  of  shock  sensitivity  (refers  to  the  ease  with  which  an 
explosive  detonates  in  response  to  die  entering  of  shock  wave),  ideally  a  large  number  of 
trails  must  be  fired  using  gap  thicknesses  near  die  estimated  critical  gap  value.  The  result  of 
such  testing,  having  the  'detonation,  non-detonation'  characteristics,  can  be  statistically 
analysed.  An  up-down  method  is  applied,  whereby,  in  a  sequence  of  trials,  the  thickness  of 
the  inert  material  or  barrier  is  increased  or  decreased  depending  on  whether  the  previous  trial 
results  in  detonation  or  non-detonation.  If  the  shot  results  in  detonation  the  next  shot  is  fired 
with  an  additional  fixed  gap  (d)  added  to  the  existing  gap  thickness,  if  it  foils  to  detonate  the 
next  shot  is  fired  with  a  removal  of  die  fixed  gap  thickness.  Usually,  twenty  to  thirty  shots 
are  required  to  complete  a  SSGT.  Using  die  Dixon  and  Mood  method  [1],  the  50%  firing  gap 
thickness,  usually  referred  to  as  critical  gap,  the  standard  error  and  the  95%  confidence  limit 
can  be  determined. 

The  standard  experimental  set-up  for  the  gap  test  [7]  at  MRL  is  shown  in  Appendix 
V,  Figure  3.  It  consists  of  a  donor  charge,  an  inert  material  barrier  or  shock  attenuator,  an 
acceptor  charge  and  a  steel  witness  plate  for  detecting  detonation;  a  deep  dent  on  the  plate 
indicates  a  detonation  otherwise  a  non-detonation  is  resulted. 

4.2.2  Application  of  the  software  package. 

A  typical  test  schedule  for  a  MRL  gap  test  (SSGT)  is  shown  in  Appendix  IV,  Figure 
3;  this  includes  the  results  for  the  statistical  analysis  of  the  data,  which  is  comparable  with  the 
results  produced  by  the  hardcopy  of  the  software  package. 

The  data  presented  in  Appendix  IV,  Figure  4  is  obtained  from  a  SSGT[7],  The 
explosive  under  test  is  PE4  with  88%  RDX  composition.  Its  density  value  is  1.59  Mg/m3 . 
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The  test  is  performed  under  ambient  temperature  using  2  mm  brass  sheet  to  make  up  the  gap 
thickness.  The  Bruceton  method  is  employed  to  determine  the  critical  gap  thickness  (m) 
require  to  detonate  the  acceptor  charge,  its  standard  error  ( Sm)  and  the  95%  confidence  limit 
(195%(m)).  The  layout  format  and  bruceton  computation  is  similar  to  that  described  in  the 
Jet  Sensitivity  Test  (Section  4. 1 .2).  The  only  difference  is  the  terminology  used  to  describe 
the  test  condition..  This  is  to  ensure  that  terminology  used  is  appropriate  for  this  application. 

4.3  Slapper  Detonator  Test 

43.1  General  description  of  test 

The  slapper  detonator  has  a  short  response  time,  is  able  to  be  precisely  timed  and 
employs  insensitive  high  explosive.  This  detonator  is  therefore  safer  to  handle  and  provides 
potentially  more  reliability  and  flexibility  than  any  other  standard  detonator. 

The  slapper  detonator  consists  of  a  capacitor,  a  flat  transmission  line,  a  fast  acting 
switch  and  a  thin  metallic  bridge  foil,  in  intimate  contact  with  the  bridge  foil  is  a  thin  sheet 
of  insulating  material,  usually  Kapton  or  Mylar,  known  as  the  flyer  plate.  Underneath  foe 
bridge  foil  is  a  tamper.  The  flyer  plate  is  separated  from  the  acceptor  explosive  by  a  short 
barrel.  A  typical  laboratory  set-up  is  shown  in  Appendix  V,  Figure  4. 

When  the  capacitor  is  discharged  rapidly,  the  necked  down  region  of  foe  bridge  foil 
vaporises.  The  flyer  plate  is  accelerated  by  foe  expanding  gas  to  a  velocity  of  several  mm  per 
ps  before  impacting  foe  acceptor  explosive.  At  impact  a  high  strength  shock  wave  is 
transmitted  into  foe  explosive.  This  shock  wave  rapidly  heats  the  explosive  near  the  point  of 
impact  which  results  in  foe  formation  of  a  detonation  wave. 

The  performance  of  the  slapper  detonator  is  determined  by  locating  foe  mean  firing 
energy,  (m)  that  produces  a  detonation  in  the  acceptor  explosive.  An  up-down  method 
[8]  is  employed  to  provide  a  E  5^  point  and  standard  deviation.  The  usual  step  size  is 
between  5  and  lOmJ.  Sample  sizes  of  about  20  shots  are  a  practical  limitation  and  are  known 
to  give  a  poor  estimate  of  foe  standard  deviation  [1]. 

The  assessment  of  the  performance  of  different  slapper  detonators  is  based  on  foe 
relative  Ejoy.  point  for  each  design.  In  addition,  all  assessments  are  related  to  a  standard 
design  in  order  to  investigate  the  effect  of  specific  design  changes. 

4.3.2  Application  of  the  software  package. 

Appendix  IV,  Figure  5  displays  a  total  of  eighteen  shots  are  subjected  to  the  Bruceton 
analysis,  whereby  the  mean  firing  energy  level  (m),  standard  error  (Sm)  and  the  95% 
confidence  limit  (  L95°Mm))  can  be  estimated.  Ideally,  a  large  sample  size  should  be 
subjected  to  Bruceton  analysis  to  insure  a  reasonable  result,  but  sample  sizes  of  about  twenty 
shots  are  a  practical  limitation.  The  bruceton  stepping  method  used  in  this  application  is 
different  from  the  Jet  Sensitivity  and  Gap  Sensitivity  test,  it  employs  the  step-down  method 
when  the  sample  fire.  This  is  apparent  on  the  numerical  table  (or  the  graphical 
representation),  it  indicates  that  at  the  first  firing  energy  level  (Test  Level),  185  mJ,  a  no-fire 
(NF)  is  resulted  and  foe  next  firing  energy  level  is  increased  to  190  mJ,  an  increment  (Step 
Size)  of  5  mJ.  At  shot  number  3,  the  firing  energy  level  of  195  mJ,  the  sample  fired  (F), 
hence  the  next  energy  level  is  reduced  to  190  mJ.  Refer  to  Section  4.1.2  for  layout 
descriptions  and  Section  8  for  symbol  definition. 


4.4  Cap  Sensitivity  Test 

4.4.1  General  description  of  test. 

MRL  is  involved  with  the  study  of  detonators,  caps  and  various  printers.  To  ensure 
that  they  function  properly  when  electrically-initiated  from  a  weapon,  they  must  pass 
stringent  electrical  tests  to  meet  RADHAS  (Radiation  Hazards).  The  set-up  for  the 
electrically-initiated  cap  sensitivity  test  consists  of  a  firing  box  and  a  firing  chamber  where 
the  cap  is  initiated  safely. 

A  typical  conducting  composition  primer  cap  is  shown  in  Appendix  V,  Figure  5. 

The  hollow  part  of  the  assembly  consists  of  an  external  cup,  internal  cup,  insulating  envelop, 
and  contact  stud  insulating  washer.  Conducting  composition  and  primer  composition  are 
pressed  in  to  die  cap  and  disc  closed.  The  capacitor  inside  the  firing  box  is  charged  up  to  a 
specific  energy  level  (up  to  500  volts)  and  is  discharged  onto  the  cap.  With  the  rapid 
discharge  time  of  die  capacitor,  a  high  current  is  released  and  flows  down  a  finite  number  of 
conducting  paths  initiating  localised  heating  (hotspots).  The  cap  functions  when  a  hot  spot 
reacts  at  die  ignition  temperature  of  the  composition. 

Test  is  carried  out  at  different  voltage  levels  to  initiate  the  device,  then  the  Bruceton 
analysis  method  is  employed  to  determine  the  mean  and  the  standard  deviation. 

Once  the  mean  is  determined,  the  first  subgroup  of  the  device  can  be  initiated  at  a  voltage 
around  this  mean.  The  process  is  repeated  for  a  number  of  voltage  levels.  The  number  of 
caps  that  fires  or  no-fires  at  each  voltage  level  is  recorded  for  Probit  or  AMCR  analysis.  This 
process  gives  an  estimate  of  die  normal  distribution  performance  of  the  device. 

4.4J2  Application  of  the  software  package. 

Appendix  IV,  Figure  6  depicts  a  probit  printout.  A  series  of  five  voltage  levels 
incremented  by  fifty  volts  (the  increments  do  not  have  to  be  constant)  is  used.  A  total  of 
twenty  shots  are  used  at  each  level  with  the  number  of  fires  and  no-fires  recorded;  this  is  also 
represented  graphically.  The  Empirical  Probit,  Provisional  Probit  and  Expected  Probit  are 
calculated  for  cadi  voltage  level.  The  "Statistical  Results"  box  lists  the  percentages  (%) 
corresponding  to  the  mean  sensitivities  (Mx)  and  die  standard  errors  (Sx)  which  are  in  non¬ 
log  or  log  form  depending  on  whether  the  original  or  the  log-transformed  data  are  used.  The 
95%  left-hand  confidence  limits  (x95%)  are  also  listed.  The  "Mean  &  StdDev.  Results  for  x" 
box  displays  the  mean  (m),  the  standard  deviation  of  the  mean  (Sm),  the  standard  deviation 
(sd)  and  the  standard  deviation  of  sd  (Ss).  The  results  in  this  box  are  relevant  to  the  values  in 
the  "x"  column  on  the  numerical  table  (i.e.  log  or  non-log). 

Appendix  IV,  Figure  7  depicts  an  AMCR  printout.  The  caps  are  subjected  to  eight 
voltage  levels  with  a  fifty  wilt  increment.  The  AMCR  method  requires  a  constant  increment. 
The  first  level  is  the  first  level  at  which  all  the  samples  failed  to  fire  and  the  last  level  is  the 
level  at  which  all  the  samples  fired.  The  number  of  fires  and  no-fires  are  listed  together  with 
the  decimal  fractions  of  the  misfires  (Pi),  the  variance  factors  (Ki)  and  the  KiPi  values.  The 
"Statistical  Results"  box  contains  the  mean  (m),  the  standard  deviation  corrected  for  test 
interval  (sd),  the  standard  deviation  uncorrected  for  test  interval  (sd*).  The  two  other 
parameters  listed  are  the  values  of  KiPi  and  Pi  summed  over  die  levels.  Each  column  on  die 
bar  chart  represents  the  proportion  of  fires  at  each  level. 


5.  Conclusion 


A  computer  software  package  has  been  developed  for  the  statistical  analysis  of  data 
from  the  testing  of  explosive  materials  and  devices.  The  package  includes  the  choice  of  three 
statistical  methods;  Bruceton,  Probit  and  AMCR.  Advantages  of  the  scheme  include  the 
rapid  and  accurate  processing  of  data,  reduction  in  human  errors,  and  the  printout  is 
presentable  and  informative. 

The  package  has  been  successfully  adapted  to  and  trialed  with  four  small  scale 
explosive  tests. 
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8.  Glossary  of  Terms 
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is  a  pop-up  window  indicating  that  the  application  expects  some  input 
from  the  user.  It  is  closely  related  to  displaying  a  window. 

is  a  list  of  application  items  or  commands  that  the  user  can  select  to 
perform  an  operation. 

is  a  window  filled  with  a  list  of  items  that  the  user  can  select  from. 

is  a  control  button  (usually  circular  in  shape)  that  offers  the  user  a 
selection  and  represents  the  state  of  the  selection;  ON  or  OFF.  Radio 
buttons  are  usually  grouped  to  provide  die  user  with  a  choice  of  mutually 
exclusive  buttons. 

is  an  area  assigned  to  enable  user  to  input  text  to  be  used  by  die  program. 

is  a  control  indicating  the  relative  position  of  a  specific  range  along  the 
List  Box  or  Window.  It  can  be  either  horizontal  or  vertical. 

mean  of  the  sample  tested. 

standard  deviation  of  the  sample  tested. 

sample  standard  deviation/error  of  the  mean. 

mean  of  the  population. 

population  standard  deviation. 

is  given  by  the  t  distribution  table  for  sum  of  sample  tested. 

test  levels,  ie.  cover  thickness,  voltage  level,  gap  thickness  etc. 

represents  the  Empirical,  Provisional  and  Expected  Probit. 

observed  percentage 

mean  of  a  given  percentage  threshold. 

standard  deviation  of  Mx  or  Log  Mx 

is  the  95%  confidence  limit,  ie.  ±t0  95  <sm. 

95  percent  confidence  limit  of  the  mean,  ie.  m±t095  am 
95  percent  no-fire  confidence  limit,  ie.  x  %  - t095  am. 
standard  deviation  not  corrected  for  test  interval, 
refer  to  equation  2. 
refer  to  equation  3. 
refer  to  section  2.1,  C 
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Appendix  I :  Program  Structured  Diagrams 


Figure  1.  Main  Structure  Diagram 


Figure  2.  File  Menu  Structure  Diagram 
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Appendix  II :  Bruceton  Reference  Graphs 

Graphs  used  in  Dixon  and  Mood  [1] 
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Appendix  111 :  Input  and  Output  Data  File  Format 


All  input  and  output  data  files  are  saved  as  standard  ASCII  format,  so  the  content  of  the  file 
can  be  viewed  and  modified  outside  the  software,  using  other  editor  packages  such  as  Turbo 
Pascal,  qEdit,  XTGold,  etc. 

Bruceton  Data  File  Format  : 


Line  1 


Line  2-15 


Line  16 


Line  17 
Line  18...n 


Linen+1 

Linen+2... 


represent  the  file  identifier  and  condition  flags  for  the  computer 
model. 

-  variable  1  is  the  statistical  data  file  identifier 

-  variable  2  is  the  method  identifier 

-  variable  3  is  the  experimental  detail  format  identifier 

-  variable  4  is  the  terminology  identifier 

-  variable  5  is  the  system  unit  identifier 

-  variable  6-1 1  is  not  used.  It  is  there  for  future 
enhancements  without  the  need  to  convert  old  data  format. 

represent  the  environment  and  condition  in  which  the  tests  are 
conducted.  These  details  correspond  to  field  names  and  information 
in  the  computer  model  (Section  3.3.3). 

"INC"  or  "DEC"  is  the  flag  to  determine  the  type  of  stepping 
methods. 

"INC"  represents  an  increment  by  a  step  interval  on  the  test  level, 
when  a  'fire'  is  resulted,  while  "DEC"  decrements  the  test 
level  by  a  constant  step  interval  when  a  ’fire'  is  detected, 
represent  die  number  of  Bruceton  shots,  step  interval  value  and  the 
number  of  ranging  shots  respectively, 
represent  the  Bruceton  series. 

-  variable  1  is  die  test  level, 

-  variable  2  is  the  result;  T  is  'Det.'  and  'O'  is  'NonDet'. 

-  variable  3  is  the  shot  cross  reference  detail, 
represent  the  first  ranging  shot  in  a  series 
represent  the  second,'  third,  etc. 


Prabit  and  AMCR  Data  File  Format . 


Line  1  As  in  Bruceton  format. 

Line  1-15  As  in  Bruceton  format. 

Line  16  number  of  test  levels  in  the  series. 

Line  16+n  represents  the  series  of  probit  data. 

-  variable  1  is  the  test  level 

-  variable  2  is  the  number  shot  results  as  a  'fire' 

-  variable  3  is  the  number  of  shot  results  as  a  'no-fire* 
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Appendix  IV:  Application  Results  and  Outputs 
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Figure  1:  A  typical  test  schedule  for  MRL  Jet  Sensitivity  Test. 


DST04* 

AU  «TWALIA 


MRL  -  BRUCBTON  ANALYSIS 


Title  :  Covered  Composition  B 

Oate  :  8/80-1/82  OIC  :  Learmonth/Chick 


Type  :  Composition  8 

Composition  :  RDX/TNT/Wax  58/45/1 

Fabrication  :  Cast  into  cylinders  38mm  dia.  x  50mm 


Density  :  1.65  g/cm3 
Addjnfo  :  Charge  Length  100mm 


1 1 fT  I  1 1 « W  IT  3 


Set  Up  :  C 

Type  of  Shape  Charge  :  MRL  standard  38mm 
Barrier  Material  :  Mild  Steel 
Add  Info  : 


Shot  Details 


SCS21 
SC567 
SC522 
SC5 
SC524 
SC589 
SCS26 
SC570 
SC527 
SC576 
SCS28 
SC577 
SC522 
SC578 
SC529 
SC579 


SC580 


Test  Level  Result  No 


63.600 

67.200 

63.600 

57.200 

63.600 

67.200 

63.600 

57.200 

63.600 

57.200 

63.600 

57.200 
50.800 

57.200 

63.600 

57.200 

63.600 

57.200 


No. 

Shot  Details 

Test  Level 

Result 

19 

SCS86 

63.600 

H 

20 

SC584 

57.200 

21 

SC593 

63.600 

ND 

Level 

63.800 


Step  Size  :  6.400 


e  o«t. 


•  NonOot  I  : 


10  12 
Shot  No. 


Mean.m  :  59.760 
Sm  :  0.839 
L95%(m)  :  61.886.  57.634 


M  :  0.090 
•  :  0.298 
*d  :  1.907 


G  :  1.392 
195%  :  2.128 


Figure  2a:  Jet  Sensitivity  Test  on  covered  Comp  B. 
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DST04* 

AUgTgjALjA 


M  R  L  -  BRUCETON  ANALYSIS 


Title  :  Covered  H-6 

Date  :  OCT-OEC  85  OiC  :  Wiese/Bussell/Chick 


Type  :  H-6 

Composition  :  RDX/TNT/AI/Wax  45/30/20/5 
Fabrication  :  Cast  into  cylinders 

Density  :  1.739  g/cm3 

Add  Info  : 


Set_Up  :  C 

Type  of  Shape  Charge  :  MRL  standard  38mm 
Barrier  Material  :  Mild  Steel 
Add  Info  : 


Shot  Details  Test  Level  Result 


Mean.m  :  72.143 

Sm  :  2.619 

L95%(m)  :  79.650,  64.635 


Figure  2b:  Jet  Sensitivity  Test  on  covered  H-6. 


M  :  0.894 
(  :  1.171 
sd  :  7.026 


G  :  0.986 
L95%  :  7.607 
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typical  test  schedule  for  MRL  Shock  Sensitivity  Test. 


DST04t 

AUSTWAUA 


MRL  -  BRUCETON  ANALYSIS 


Title  :  Gap  Test  Assessment  of  Shock  Sensitivity 

Date:  12-2-82  OIC  : 


TYPE  :  PE4 

COMPOSITION  :  88%  ROX 

DETAILS  OF 

FABRICATION 

OENSITY  :  1.59  Mg/m  *3 

ADD.  DATA  :  Lot.  No.  80,  Incorp.  No.  690 


TEST  TEMPERARURE  :  AMB 
TYPE  OF  DONOR  :  Scale  1 

BARRIER  DETAILS  :  Brass,  2mm  thickness  interval 
Add.  INFORMATION  : 


Test  Level  Result  No 


20.000 

22.000 

24.000 

22.000 

20.000 

22.000 

20.000 

22.000 

20.000 


No. 

Shot  Details 

1 

1 

Shot  Details 


10 

10 

22.000 

NO 

11 

11 

20.000 

0 

12 

12 

22.000 

NO 

13 

13 

20.000 

0 

14 

14 

22.000 

0 

15 

15 

24.000 

NO 

16 

16 

22.000 

NO 

17 

17 

20.000 

NO 

18 

18 

18.000 

0 

Mean,m  : 
Sm  : 
L95%(m) 


Test  Level 


20.000 

22.000 

24.000 

22.000 

24.000 

22.000 

24.000 


Result 


21.667 

0.419 

:  22.681,20.652 


M  :  0.389 
s  :  0.677 
sd  :  1.354 


G  :  1.072 
L9S%  :  1.014 


Figure  4:  Shock  Sensitivity  Test  on  PE-4,  using  customised  title  format. 
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DST04* 


M  RL  -  BRUCETON  ANALYSIS 


Title  :  HNS  42 

Data  :  1990 


TYPE  :  HNS  42 
SSA  :  "42  m*2/g 
DENSITY  :  90%  TMD 
ADD.  DATA  : 


SET-UP  :  Standard  Slapper  Detonator 
TYPE  OF  DEVICE  USED  :  0. 1 89uF  capacitor 
ADD.  INFORMATION  : 


No. 

Shot  Details 

1 

1 

Test  Laval  Result  No. 


185.000 

190.000 

19S.000 

190.000 

185.000 

190.000 

185.000 

190.000 

195.000 

200.000 

195.000 

190.000 

195.000 

190.000 

185.000 

180.000 

185.000 

180.000 


Shot  Details  Test  Laval  Result 


Laval 

200.000 

195.000 

190.000 

185.000 

180.000 


Step  Size  :  5.000 


•  Rro 


lloft,  I  FWtanw  :  figS.bru 


8  10 
Shot  No. 


Mean.m  :  190.000  M  : 

Sm  :  2.171  •  : 

L95%(m)  :  195.904,  184.096  -  : 


Figure  5:  Slapper  Detonator  Test,  using  customised  title  format. 


G  :  0.973 
195%  :  5.904 
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DST04* 

AUgTWAUA 


M  R  L  -  PROBIT  ANALYSIS 


Title  :  n43  Primers  Stressed 

Oate  :  8/6/90  OIC  :  L.  Redman 


TYPE  :  CAPS 
COMPOSITION  :  CC 
ADD.  DATA  : 


SET-UP  : 

TYPE  OF  DEVICE  USED  : 
ADD.  INFORMATION  : 


100.000 

150.000 

200.000 

260.000 

300.000 


%  Fire 


5.000 

20.000 

40.000 

45.000 

75.000 


Exp.Probit 


Statistical  Results  : 


Mean  &  StdDev.  Results  for  x 

Mean.m 

:  2.370 

Sm 

:  0.033 

sd 

:  0.222 

S* 

:  0.049 

DST04* 

AUTWAtIA 


MRL  -  AMCR  ANALYSIS 


Tit!a  :  N43  Caps  in  Primers  Unstressed 

Oats:  20/11/91  OtC  :  L.Redman 


TYPE  :  CAPS 
COMPOSITION  :  CC 
ADO.  DATA  : 


SET-UP  : 

TYPE  OF  DEVICE  USED  :  Caps 
ADD.  INFORMATION  : 


Laval 

fire 

100.000 

0 

160.000 

1 

200.000 

3 

260.000 

10 

300.000 

14 

360.000 

18 

400.000 

18 

450.000 

20 

Statistical  Results  :  | 

Sum(KiPi) 

:  18.700 

Sum(PI) 

:  3.800 

Mean.m 

:  315.000 

sd 

:  75.166 

sd* 

:  1.503 

FileName  :  fig7.amc 


Figure  7:  Conducting  Composition  Cap  Sensitivity  Test,  using  customised  title  format. 


Appendix  V:  Experimental  Test  Set-up 
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V>. 


M  R  L  SCALE  I  DONOR 


Paper  disc 


34  mm 


Pressed  p6TN 


12.7  mm 


BRASS  LAMjNATg 


wiouioea  pusuc  r»«*o  * - . - 

Bridge  wire  ^ 

-i 

-*!  H 

••  i* 

Hi — w-i 

jh 

**  ** 
li 

Low  density  _ 

\  \ 
x\ 

_ _ [} 

precipitated  PETN 

0*. 

' 

Pressed 

33 

MRL  SCALE  2  DONOR  bigh  density  PE  TN 

*  •:  . 

'Vs  . 

Donor 


Laminated 
brass  barrier 
(gap) 


Acceptor 

explosive 


Mild  steel 
witness  block 


Figure  3:  The  MRL  Small  Scale  Gap  Test  (SSGT)  assembly. 
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Figure  4:  Slapper  Detonator  Assembly. 


Figure  5:  A  typical  Conducting  Composition  Cap. 
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